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Abstract 

The experimental Raman spectra of graphene exhibit a few intense two-phonon bands, which are 
enhanced through double-resonant scattering processes. Though there are many theoretical papers 
on this topic, none of them predicts the spectra within a single model. Here, we present results 
for the two-phonon Raman spectra of graphene calculated by means of the quantum perturbation 
theory. The electron and phonon dispersions, electronic lifetime, electron-photon and electron- 
phonon matrix elements, are all obtained within a density- functional-theory-based non-orthogonal 
tight-binding model. We study systematically the overtone and combination two-phonon Raman 
bands, and, in particular, the energy and polarization dependence of their Raman shift and inten- 
sity. We find that the ratio of the integrated intensities for parallel and cross polarized light for all 
two-phonon bands is between 0.33 and 0.42. Our results are in good agreement with the available 
experimental data. 
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I. INTRODUCTION 



Presently, graphene is considered as a prospective material for nanoelectronics and 
nanophotonics.-*^ Among the various experimental techniques, the Raman spectroscopy has 
proven to be an indispensable tool for investigation of this material.-"- 

Graphene has a single Raman-active phonon observed as an intense line (the G 
band) in the first-order Raman spectra. The second-order spectra of graphene with low 
defect density has several intense bands, which originate from scattering of electrons and 
holes by two phonons of the same/different frequency and non-zero momentum and are 
called overtone/combination bands. The appearance of intense second-order bands can 
be explained by the double-resonant (DR) scattering mechanism.-i^ These bands contain 
valuable information on the phonon dispersion^i"— and the electron-phonon and electron- 
electron matrix elements.™ The Raman spectra of graphene with defects show additional 
bands, which arise from DR scattering of electrons and holes by phonons and defects. 

The theoretical investigation of the two-phonon DR scattering in graphene has been 
performed using various approximations: replacement of the electron-photon and electron- 
phonon interactions with constants, using a constant electronic lifetime, considering only 
high-symmetry directions in the Brillouin zone, as well as exact DR conditions. The pre- 
dicted dispersive behavior of the Raman bands'' and the frequency shift of the Stokes and 
anti-Stokes Raman bands^^ have been found in quantitative agreement with the experimental 
data. It has also been realized that the integration over the entire Brillouin zone of graphene 
is essential for predicting the Raman intensity.— While most of the theoretical papers fo- 
cus on the most intense overtone band, in a recent study, all two-phonon Raman bands 
with observable intensity have been calculated using an electron-phonon matrix element de- 
rived within a nearest-neighbor vr-band tight-binding model.— The dominant contribution 
to the two-phonon bands from different parts of the Brillouin zonei2.ii^ and from different 
scattering processes-- has also been discussed. The polarization dependence of the most 
intense overtone band has been studied experimentally and theoretically.— The progress, 
made so far in the modeling of the DR bands, has been achieved either with simple vr-band 
tight-binding models, or with more sophisticated models but relying on approximations of 
the electron-phonon matrix element, and in most cases concerns only a few intense bands. 
Our experience indicates that, although the Raman shift of the bands is not sensitive to the 
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used matrix element, their Raman intensity crucially depends on it. 

Here, we calculate the two-phonon DR Raman spectra of graphene using a non-orthogonal 
tight-binding (NTB) model, which implements parameters derived from a density functional 
theory (DFT) study and thus has no adjustable parameters.— In particular, the electronic^ 
and phonon^^ dispersion, the electron-photon and electron-phonon matrix elements,— as 
well as the electronic linewidth^i are all obtained within this model. The NTB model for 
electrons and phonons in graphene is introduced in Sec. II. The calculated two-phonon DR 
Raman spectra of graphene and its polarization dependence are discussed in Sec. III. The 
paper ends up with conclusions (Sec. IV). 

II. THEORETICAL PART 
A. The NTB model 

We use a NTB model with four valence electrons per carbon atom to calculate the elec- 
tronic dispersion of graphene. This model is based on matrix elements of the Hamiltonian 
and overlap matrix elements derived from DFT— and therefore it does not rely on any ad- 
justable parameters. It also allows one to estimate the total energy and the forces on the 
atoms. This feature is utilized for relaxation of the atomic structure. Up to a few electron 
volts away from the Fermi energy, the electronic structure of graphene has the form of conic 
valence and conduction bands (Dirac cones) with a common apex (the Dirac point) at two 
non-equivalent special points, K and K', of the Brillouin zone. This specific form of the 
electronic bands plays an important role in the enhancement of the two-phonon Raman 
scattering through the DR mechanism. 

The dynamical model of graphene uses a djTiamical matrix derived by a perturbative 
approach within the NTB model.— The electron-photon and electron-phonon matrix ele- 
ments are calculated explicitly.— The summation over the Brillouin zone in the first-order 
perturbation term of the dynamical matrix is performed over a 40 x 40 mesh of k points, 
for which the phonon frequencies converge within 1 cm~^. The calculated in-plane phonon 
branches of graphene, after scaling by a factor of 0.9, agree fairly well with the available 
experimental data^^ (Fig- !)• The phonons with displacement in the graphene plane (in- 
plane phonons) can interact with electrons and thus can contribute to the Raman spectra. 
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Those with atomic displacement perpendicular to graphene (out-of-plane phonons) are less 
well reproduced but they do not contribute to the spectra. 

It will be shown below that only phonons, close to the high-symmetry directions FK, FM, 
and KM of the Brillouin zone, are of major importance for the two-phonon spectra. The 
phonon branches along these directions will be denoted, as usual, by two-letter acronyms 
describing their vibrational pattern: the letters O and A stand for "optical" and "acoustic" , 
respectively; the letters L, T, and Z denote in-plane longitudinal, in-plane transverse, and 
out-of-plane atomic displacement, respectively. The acronyms for the branches along the KM 
direction will be primed. Alternatively, for each wavevector, the phonons with be ascribed 
the index v = 1,...,6, in order of increasing frequency. The phonons with a certain u 
can belong to branches with a different vibrational pattern. For example, a phonon with 
z/ = 6 can belong to the LO or TO branch. It will also be argued that only phonons close to 
the F and K points give a significant contribution to the two-phonon Raman spectra. Such 
phonons will be denoted by acronyms ending with @F and @K. For example, LO phonons 
close to the F point will be denoted by LO@F and TO phonons close to the K point along 
the KM direction will be denoted by TO'@K (see. Fig. 1). 



B. The double-resonant processes 

The amplitude for two-phonon DR Raman scattering processes in graphene is described 
by fourth-order terms in perturbation theory.^- The underlying processes include virtual 
scattering of electrons/holes by phonons between states of the Dirac cone at the K point or 
the K' point, or between states of the Dirac cones at the K and K' points. The momentum 
is conserved in each virtual process but the energy is conserved only for the entire DR 
process. Below, we will consider only Stokes processes. In this two-phonon DR 

process includes an absorption of a photon with a creation of an electron-hole pair, two 
consecutive processes of scattering of an electron/hole with creation of a phonon, and a 
recombination of the electron- hole pair with an emission of a photon (Fig. 2). There are 
altogether eight such processes.— The total two-phonon Raman intensity is given by the 
expression 
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Here, the inner sum is the scattering amphtude. Eu, u = i,a, b, c, /, are the energies of 
the initial [i), intermediate (a,6, c), and final (/) states of the system of photons, electrons, 
holes, and phonons. In the initial state, only an incident photon is present and, therefore, 
Ei = El, where E^ is the incident photon energy. In the final state, there is a scattered 
photon and two created phonons. M^v are the matrix elements for virtual processes between 
initial, intermediate, and final states. In particular, Mai and Mj^ are the matrix elements 
of momentum for the processes of creation and recombination of an electron-hole pair, re- 
spectively. Mf,a and Met are the electron/hole-phonon matrix elements. 7 is the sum of the 
halfwidths of pairs of electronic and hole states, and will be referred to as the electronic 
linewidth. The electron-photon and electron-phonon matrix elements, and the electronic 
linewidth are calculated explicitly.— The Dirac delta function ensures energy conservation 
for the entire process. In the calculations, it is replaced by a Lorentzian with a halfwidth of 
5 cm~^. The summation over the intermediate states runs over all valence and conduction 
bands, and over all electron wavevectors k. The summation over the final states runs over all 
phonon branches and phonon wavevectors q. For both summations, convergence is reached 
with a 800 x 800 mesh of k and q points in the Brillouin zone. 

For the discussion of the polarization dependence of the Raman intensity of the two- 
phonon bands, it is advantageous to rewrite Eq. ([1]) in the form 



Here, e^, and 65 are the polarization vectors of the incident and scattered laser light, re- 
spectively, and R is the Raman tensor. Everywhere below we consider only backscattering 
geometry in accord with the usual experimental Raman setup for graphene and, therefore, 
the polarization vectors lie in the graphene plane. 
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III. RESULTS AND DISCUSSION 



A. Electronic linewidth 

The electronic linewidth 7 is due to a large extent to scattering of the electrons (holes) by 
phonons and other electrons (holes). The majority of the published reports assume that 7 is 
energy-independent. Recently, it has been argued that in undoped graphene 7 is dominated 
by electron-phonon processes and the expression 7 = 9A4E + SAOE^ has been derived, 
where the energy separation between the valence and conduction bands E is in eV and 7 
is in meV.-i^ Here, the electronic linewidth was calculated by summing up the contributions 
of all electron/hole-phonon scattering processes for all phonons in the Brillouin zone as a 
function of E. The obtained energy dependence was approximated in the range [1.0, 3.5] eV 
with the expression 

7 = 12.60E + 3.45^2 (3) 

For energies in this energy range, 7 changes more than four times from 16 to 86 meV (Fig. 
3). 

B. Overtone bands 

The calculated overtone Raman spectrum for El = 2.0 eV and parallel light polarization 
along a zigzag line of carbon bonds is shown in Fig. 4 (bottom). It has two intense bands 
(2D and 2D') and three weaker ones {2D^, 2D^, and 2D") and a much weaker one (2D^). 
There are also other bands of in-plane phonons close to the F and K points but they are 
either very weak, or are in the shoulders of intense bands, and in both cases are practically 
unobservable. 

The assignment of the overtone spectrum can be performed by analyzing the contribution 
of phonons with different u and from different parts of the Brillouin zone of graphene. First, 
the contributions of the phonons with u = 2, ...,6 are given in Fig. 4 (top five graphs). The 
bands in these spectra originate from pairs of phonons TA@r [u = 2), LA@r and TA@K 
(z/ = 3), LA@K (z/ = 4), TOOK (z/ = 5), and TO'@K and LO@K {u = 6) (see also Fig. 
1). The contributions for u = 5,6 are by about three orders of magnitude larger than the 
remaining ones. 
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Secondly, the assignment of the Raman bands to definite phonons at the T and K points 
is supported by the analysis of the contributions to the spectra from different parts of 
the Brillouin zone for parallel light polarization with averaging over all orientations in the 
graphene plane. In Fig. 5, the regions with major contribution to the bands are given 
by shaded areas. In particular, the 2D^ band comes from phonons TA@r along the FK 
direction, the 2D'^ band - from phonons LA@F along FM, the 2D" band - from phonons 
LA@K along FK (not shown in Fig. 5), and the 2D' band - from phonons LO@F from all 
directions in the Brillouin zone. The 2D band originates from phonons TO@K and TO'@K, 
which contribute to the integrated Raman intensity of this band in ratio ^ 10 : 1. It has 
long been accepted that the dominant contribution to the 2D band comes from phonons 
along the KM direction and that the contribution along the FK direction vanishes because 
of interference.— Our graphs in Fig. 4 and 5 for u = 5,6 confirm the recent conclusion that 
the main contribution to this band comes from the FK direction and only a small part comes 
from the KM direct ion. ■ii'=i^ We also confirm the result of Ref.— , that among the eight DR 
processes, those with scattering of an electron and a hole have dominant contribution to the 
bands. 

The overtone bands are dispersive, i.e., their Raman shift depends on the laser energy 
(Fig. 6). The shift increases (decreases) due to increase (decrease) of the phonon frequency 
away from the F and K points. It is approximately linear in El with slopes of ~ for the 
2D' band, 89 cm-^/eV {2D band), -193 cm'VeV {2D"), 419 cm-7eV {2D^), and 256 
cm~^/eV (2D^) at = 2.0 eV . The slope of the 2D band corresponds to the experimental 
one of 88 cm-VeV (Ref. fii]). 

The integrated Raman intensity of the overtone bands (Fig. 7) is quasi-linear in El 
except for the 2D^ band. The ratio of the integrated intensities A {2D) /A (2D'j has been 
discussed a lot in the literature because it is related to the electron-phonon matrix elements 
at the K and F points. Here, we find that this ratio depends on the energy: it varies from 
4.7 for El = 1.0 eV to 15.0 for El = 3.5 eV; for El = 2.4 eV it is equal to 8.7. By contrast, 
recent sophisticated DFT-GW calculations!^ yield 21.5, which is several times larger than 
our value. The reason for this disagreement can be found in the use of GW corrected 
electron-phonon matrix elements for the LO phonon at the F point and the TO phonon 
at the K point. While the ratio of the squares of these matrix elements at the two points 
derived by DFT is M^/M^ = 2.02, the DFT-GW result is 3.03 (Ref. f^]). Since the Raman 
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intensity depends on the square of this ratio, the DFT-GW gives an intensity ratio that is 
larger than the DFT one by a factor of 2.25. Our resuh, corrected by the same correction 
factor, is 19.6. Both theoretical results underestimate the experimental values of 26 (Ref. 
P^]) and 27 (Ref. [^]) by ~ 25%. The origin of this underestimation is still unknown. 

C. Combination bands 

The combination Raman spectrum is calculated for El = 2.0 eV and parallel light polar- 
ization along a zigzag line of carbon bonds (Fig. 8, bottom). The spectrum is dominated by 
two intense bands: D + D" and a higher frequency one, which originate from the branches 
5 + 4 and 6 + 5, respectively. There are also four other less intense bands: D' + D^, D' + D^, 
and D + D'\ 

As above, the assignment is facilitated by considering the contributions to these bands 
from pairs of phonons with i/, z/' = 2, 3, 4, 5, 6. It is seen in Fig. 8 (top six graphs) that for 
each pair v + u' there is a single band. These bands are due to pairs of phonons TOTA@K 
{v + v = 5 + 3), TOLAOK (5 + 4), LOTAOF (6 + 2), LOLAOF (6 + 3), T0'L0'@K (6 + 4), 
and TO LA @K (6 + 5) (see also Fig. 1). Unlike the case of the overtone bands, the intensity 
of the various combination bands varies only by one order of magnitude. The band D + D" 
is mainly due to phonons TOLA@K with a small contribution of phonons T0'L0'@K. The 
phonons TO LA @K give rise to an intense band at ~ 2700 cm . The large contribution of 
the phonons T0'L0'@K and T0'LA'@K is consistent with the large electron-phonon matrix 
element for the phonon TO'@K and the nonzero matrix element for the phonons LA'@K and 
LO'@K. We note that the latter contrubite to the overtone spectra as well but their bands 
are masked in the shoulder of the intense 2D band. The band at 2700 cm~^ also overlaps 
considerably with the much more intense 2D band and cannot be observed as a separate 
feature. Thus, the total Raman spectrum exhibits the intense 2D band, the two weak bands 
2D' and D + D" {Fig. 9), as well as the two very weak bands D' + D^ and D' + D^. 

Similarly to the overtone bands, the combination bands are dispersive. The dependence of 
the Raman shift on the laser energy (Fig. 10) is almost linear with a slope of 74 cm~^/eV for 
the highest-frequency band, —50 cm~^/eV {D + D"), —54 cm~^/eV {D + D^), 210 cm~^/eV 
{D' + -D^), and 123 cm~^/eV {D' + D^) at El = 2.0 eV. Each slope is equal approximately 
to the sum of half of the slopes of the corresponding overtone bands. The calculated energy 
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dependence of the D + D band frequency reproduces fairly well the experimental data but 
the derived slope is almost three times larger than the experimental one of —18 cm~^/eV 
for energies between 1.92 and 2.71 eV (Ref. ^])- The calculated frequencies for the D' + 
and D' + bands underestimate the Raman data up to 80 and 30 cm^^, respectively. 
Similar underestimation is present in other precise calculations.—"^ On the other hand, the 
frequency slopes of the bands D' + and D' + are in fair agreement with the measured 
ones of 221 cm~^/eV and 140 cm~^/eV, respectively.— We note that neither we, nor the 
authors of Ref. fi^], have found any observable combination band for phonons TOLA@r 
though such band has been establishect^ by fitting a low-intensity Raman band with two 
Lorentzians for two overlapping bands assigned to phonons TOLA@r and LOLA@r. 

The calculated integrated Raman intensity is shown in Fig. 11. The curves are increasing 
functions of laser energy for all bands except for the band D' + and the highest frequency 
one. For the most studied D + D" band, we find A {2D) /A {jD + D"^ = 12.2. Introducing 
the DFT-GW correction factor for the electron-phonon matrix element of the TO phonon at 
the K point of 1.76 (Ref.-^), for the latter ratio we obtain the value of 21.5, which agrees with 
the previous estimate^"^ of 16 and the experimental one'^ of ~ 21. The inset of Fig. 11 shows 
the dependence of the DFT-GW corrected ratios A {2D') /A {2D) and A (^D + D") /A {2D) 
as a function of E^. Both curves agree well with the theoretical ones^ and the available 
experimental data.- 

D. Dependence on the electronic linewidth 

Most of the previously reported two-phonon Raman spectra have been calculated for 
energy-independent 7. In order to study the effect of this approximation we performed 
calculation of the integrated intensity Ac for constant linewidth 7c = 7(-El = 2.0 eV) = 39 
meV as well. We found that, for any of the two-phonon peaks discussed above, the calculated 
ratio of the integrated intensities, A and A^, for variable and constant linewidth, respectively, 
can be fitted very well with the expression A/ Ac = (7c/7)^- Therefore, A can be written as 
A = //7^, where / depends on El and does not depend on 7. It is clear from Figs. 7 and 
11 that the rate of change of / as a function of El is different for the different two-phonon 
bands. Thus, the ratio of the integrated intensities for any pair of two-phonon bands 1 and 
2, A1/A2 = /1//2, should depend on El in agreement with the recent conclusions.— We 
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note that, in the latter paper, deviation from the 1/7^ behavior has been found. The correct 
dependence should be established by comparison with experimental Raman data collected 
at more values of the laser energy. 

The intensity ratio for the 2D and 2D' overtone bands can be obtained using the simple 
tight-binding analytical formulas:— A(2D) oc 2(77^/7)^ and A{2D'^ oc (7r/7)^. Here, 
7 = 7r + 7i^ is the total electronic linewidth; 7r,i^ oc E^M^ j^/ur^K are the linewidths due 
to the LO phonon at the F point and the TO phonon at the K point, Mp^, and ur^x are 
the corresponding square of the electron-phonon matrix element and the phonon frequency, 
respectively. The predicted ratio of the integrated intensity of the two bands is a constant, 
equal to 3 (or 7.9), evaluated with electronic linewidth from DFT (or DFT-GW).— The 
latter values are about 7 (or 2.7) times smaller than those of the precise derivations at 
El = 2.4 eV, presented here and in Ref. The disagreement for the intensity ratio can 
be sought in the approximations used in the theoretical scheme of Ref. f^]. 

E. Polarization dependence of the bands 

Let us choose a coordinate system in the graphene plane with x axis along a zigzag line 
of carbon bonds (and therefore the y axis is along an armchair line of carbon bonds), and 
set ei = (cos a, sin a) and = (cos /3, sin /3). In the case of parallel and cross polarizations, 
the intensity is not angle-dependent because of the high-symmetry of graphene. On the 
contrary, for fixed polarization angle and variable analyzer angle, or vice versa, the intensity 
depends only on the difference a — P and does not depend on the sign of this difference. 
Therefore, the Raman intensity, Eq. can be written as 

I = {l\\-I±)cos^{a- (3) + I^ (4) 

The intensity for parallel and cross polarized light, /|| and I±, can be obtained by fitting 
this expression to the calculated angular dependence of the intensity. Expression, similar to 
Eq. dl]), holds for the integrated intensity as well. 

Fig. 12 shows the polarization dependence of the ratio of the calculated integrated 
intensities y4_|_/A|| at E^ = 2.0 eV. As could be expected, all two-phonon bands have similar 
angular dependence of the Raman intensity, which follows from Eq. The derived values 
of the ratio A^/A^ are: 0.3678 {2D^), 0.4195 {2D^), 0.3321 {2D"), 0.3377 {2D), 0.4075 
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{2D'), 0.3718 (D' + D^), 0.4073 {D' + D"^), 0.3401 (D + D^), 0.3350 {D + D"). The intensity 
ratio for all bands is between 0.33 and 0.42. The ratio for the 2D band corresponds to the 
value 1/3 derived within a simple tight-binding model^ and the tight-binding estimate^ of 
0.32, and is in fair agreement with the experimental value^^ of 0.3. So far as we are aware, 
there are no reports on the polarization dependence of the other two-phonon bands. 

IV. CONCLUSIONS 

We have presented a complete theoretical treatment of the two-phonon Raman bands of 
graphene within a non-orthogonal tight-binding model with no adjustable parameters and 
no other approximations. We have calculated the laser energy and polarization dependence 
of the Raman shift and intensity. In particular, the ratio of the integrated Raman intensity 
for parallel and cross polarized light for all bands is between 0.33 and 0.42. The agreement 
of our results with available experimental data is very good and the predictions can be 
used for further comparison to experiment and for support in the assignment of the two- 
phonon Raman spectra. Furthermore, we have made definite conclusions on the dominant 
contributions to the Raman scattering amplitude. Our computations can easily be extended 
to include higher-order Raman processes. 
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FIG. 1. Calculated phonon dispersion of graphene along the high-symmetry directions in the 
Brillouin zone.'^^ The six phonon branches are marked by the acronyms LO, TO, ZO, LA, TA, 
and ZA. The letters O and A stand for "optical" and "acoustic", respectively; L, T, and Z denote 
in-plane longitudinal, in-plane transverse, and out-of-plane atomic displacement, respectively. The 
crosses mark the phonons which play major role in the DR processes. Inset: the hexagonal Brillouin 
zone of graphene with the special points F, M, K, and K'. 
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FIG. 2. Schematic representation of the DR processes in graphene. The sohd hnes are a cross- 
section of the Dirac cones at the K and K' points of the Brillouin zone. The dashed hnes are 
virtual processes: "eh+" and "eh~" are electron (e) - hole (h) creation and annihilation processes, 
respectively; "ep" and "hp" are electron and hole scattering processes by a phonon (p). DR 
processes can also take place between bands at the K or K' point only. 
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FIG. 3. Calculated electronic linewidth ^(E) (solid symbols). The dashed line is a fit to the 
calculated data. The fitting function is given by Eq. ([3]). 
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FIG. 4. Calculated overtone Raman spectrum for = 2.0 eV (bottom). The contributions from 
phonons with u = 2,..., 6 are shown in the top graphs. The notation of the Raman bands 2D^, 
2D^ and 2D^ is taken from Ref. 
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FIG. 5. The contribution to the overtone bands from phonons with v = 2, 3, 5, 6 from different 
parts of the rhombic Brillouin zone of graphene. The graph for = 4 for phonons 2LA@K (not 
shown) is similar to that for v = 5. 



16 



3500 
3000 
2500 
J_ 2000 
^ 1500 

CO 

E 

1000 
500 



2D' 




~0 






— • 


• 

■ 2D 
• 


— • — 
— • — 


89 ^ 


— •• 


• 

. 2D" 


■ — • — . 


-193 
■ • 








- 






• 














• 




.2D4 


• 


419 
• 

256 


• 

• — ■ 


• 


• 


' 2D3 
1 


• — 

1 


• — 

1 


1 


1 


1 



1.0 



1.5 



2.0 2.5 

El (eV) 



3.0 



3.5 



FIG. 6. Calculated dependence of the Raman shift of the overtone bands on El (solid symbols). 
The numbers are the slopes of the curves for El = 2.0 eV. The empty symbols are experimental 
data^ 




FIG. 7. Calculated dependence of the integrated Raman intensity of the overtone bands on El 
(solid symbols). 
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FIG. 8. Calculated combination Raman spectrum for El = 2.0 eV (bottom). The contributions 
from pairs of phonons with v + u' = 5 + 3, ...,6 + 5 are shown in the top graphs. 
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FIG. 9. Calculated total two-phonon Raman spectrum for El = 2.0 eV (bottom). The constituting 
combination and overtone spectra are also shown (middle and top). 
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FIG. 10. Calculated dependence of the Raman shift of the combination bands on Ej^ (solid sym- 
bols). The highest-frequency band is due to T0'LA'@K phonons. The numbers are the slopes of 
the curves for = 2.0 eV. The empty symbols - squares^^ and circleSji^fl are experimental data. 




FIG. 11. Calculated dependence of the integrated Raman intensity of the combination bands on 
El (solid symbols). The highest-intensity band at low energies is due to T0'LA'@K phonons. 
Inset: The DFT-GW corrected ratios A (2Z)') /A {2D) and A (^D + /A {2D) (small symbols) 
in comparison with experimental data^ (large symbols). The scale is the same as in Fig. 7. 
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FIG. 12. Calculated polarization dependence of the integrated Raman intensity of the 2D band 
(solid symbols). The angle a is relative the x axis along a zigzag line of carbon bonds. The 
polarizer and analyzer angles are given by the notation (a/3). The curves for parallel and cross 
polarization, [aa) and(a — 90, a), as well for fixed polarizer direction along the x axis and variable 
analyzer orientation (Oa) are shown. The solid line for polarization (Oa) is a fit of Eq. Q. The 
obtained ratio of the integrated intensities is A_l/j4|| ~ 0.34. 
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